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Abstract 

The study of laminated magnetostrictive functionally 
graded material (FGM) plate without shear deformation 
under thermal sinusoidal vibration and transient response 
is calculated by using the generalized differential 
quadrature (GDQ) method. In the thermoelastic 
stress -strain relations that containing a power-law function 
of a two-material FGM plate, the linear temperature rise 
and the magnetostrictive coupling terms with velocity 
feedback control. Four edges of rectangular laminated 
Terfenol-D FGM plate with simply supported boundary 
conditions are considered. The suitable product value of 
coil constant and control gain can be used to reduce the 
amplitude of center displacement into a smaller value. 
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Introduction 

Typical functionally graded material (FGM) is 
usually made of different phases constituent 
materials, for example, the ceramic and metal used in 
the engine combustion chamber to withstand 
ultra-high-temperature and to reduce the stress 
singularities, respectively. Chi and Chung (2006) 
presented the mechanical analysis of FGM plates 
subjected to transverse load. There were several 
studies in the transverse displacements for the 
laminated plate including the shear deformation 
effect. Amabili and Farhadi (2009) made the research 
of the shear deformable versus classical theories for 
nonlinear vibrations of rectangular isotropic and 
laminated composite plates. The over -prediction of 
natural frequencies in the solution without shear 
deformation was found. Ray and Shivakumar (2009) 
analyzed the effect of shear deformation on the 


piezoelectric fiber- reinforced composite plate by 
using the finite element method (FEM). Nguyen et al. 
(2008) obtained the static numerical results for the 
FGM plate with the effect of shear deformation. 

Magnetostrictive material Terfenol-D has the 
magneto- electric coupling property under the action 
of magnetism and mechanism. Hong (2009) used the 
computational generalized differential quadrature 
(GDQ) method to study the transient responses of 
magnetostrictive plates under thermal vibration. 
Thermal stresses and center displacement with and 
without shear effect were calculated in the thin and 
thick plate, respectively. Ramirez et al. (2006) 
obtained the free vibration solution for 
magneto-electro-elastic laminates through the Ritz 
approach. Lee and Reddy (2005) analyzed the 
non-linear response of laminated plate of 
magnetostrictive material under thermo-mechanical 
loading by using the FEM. Lee et al. (2004) obtained 
the transient vibration values of displacement for the 
Terfenol-D plate including the effect of shear 
deformation by using the FEM. Hong (2007) used the 
GDQ method to make the thermal vibration study 
for the Terfenol-D magnetostrictive laminated plate 
with the first-order shear deformation. Hong (2012) 
used the GDQ method to make the Terfenol-D FGM 
plate analyses under rapid heating induced vibration 
with the shear deformation effect. It is interesting to 
study thermal vibration of the transverse center 
displacement and thermal stress in the Terfenol-D 
FGM plate without the shear deformation effect by 
using the GDQ method. 

Formulation 

FGM 

Most materials of FGM can be used in the 
environment of higher temperature and can be 
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expressed in series form as follows by Chi and 
Chung (2006). 

Pf,.=I.P i V r 

i = 1 

(1) 

where P fgm is the material properties of FGM, n m 
is the number of materials mixed to form the FGM, 

n m 

V is the volume fractions, and > V = 1 for all 

i i 

i = 1 

constituent materials, P. is the individual 
constituent material properties, usually with the 
form as follows. 

p, = p 0 (P- + 1 + P l T + P 2 T 2 + P 3 T 3 ) . 

( 2 ) 

in which P 09 P_ X9 P X9 P 2 and P 3 are the 

temperature coefficients, T is the temperature of 
environment. 


K f gm - ( K 2 + ^1 ) / 2 • 

(3e) 

where E is the Young's modulus, v is the 
Poisson's ratio, p is the density, cc is the thermal 
expansion coefficients, /cis the thermal conductivity, 
the subscript fgm represents the FGM plate, the 
subscripts 1 and 2 represent the constituent material 
1 and 2, respectively. The property terms E x E 2 

V 1 V 2 Pi Pi a \ a 2 K \ K 1 are ex P resse d 

corresponding to term P. in equation (2). 

GDQ Method 

The GDQ method approximates the derivative of 
function, and the first-order and the second-order 
derivatives of function f*(x 9 y) at coordinates 
( x i9 y j ) of grid point (i 9 j) can be discretized in 

series forms by Shu and Du (1997) and rewritten as 
follows: 


We denote the parameters for a two-material ( n m =2) 
FGM plate as follows, a and b is the length in the 
x, y direction of the plate, h is the total thickness 
of magnetostrictive layer and FGM plate, h 3 is the 

thickness of magnetostrictive layer, h x and h 2 
are the thickness of FGM material 1 and FGM 
material 2, respectively, p x and p 2 are the in-plane 
distributed forces, q is the applied pressure load. 
The sum of volume fractions is in the form: 
V x + V 2 = 1 , the variation form of V 2 used in the 

z + hl 2 
h 

is the thickness coordinate, h is the thickness of 
FGM plate, R is the power-law index. And the 
material properties for equation (1) can be assumed 
for the simple calculation and expressed as follows 
by Hong (2012). 

(3a) 

v fgm =(v 2 +VJ/2, 

(3b) 

P fgm ~ (. Pi + A ) / 2 , 

(3c) 

a fgm =(«2 +« l )/ 2 , 

(3d) 

shear stresses and shear strains in the laminate 


power-law function is V 2 — ( 


) Rn , where z 


r 


3c 


r 


/= l 

(4a) 

M 

|. . *y B m f im , 

\i,j j ,m J l ’ m ' 

m = 1 

(4b) 
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1=1 

(4c) 
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3c 
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m = 1 


d 2 f 


2 r* N M 

Uj ~ 


dxdy 




1=1 777=1 


(4e) 

where A^ m ) and denote the weighting 

coefficients for the m th -order derivative of the 
function f*(x,y) with respect to x and y 
directions. 


Thermo Elastic Stress-Strain Relations with 
Magnetostrictive Effect 

We consider a rectangular laminated 
magnetostrictive FGM plate of the length a, b in the x, 
y direction, respectively, under uniformly 
distributed load and thermal effect. There are no 
without shear effect assumption. The plane stresses 
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Q l6 =Q 26 = 045 = 0 - 

(6d) 

P P and p are the in-plane strains, 

x * y xy 

AT = T 0 (x,y,t) + ^rT l (x,y,t) is the temperature 

h 

difference between the FGM plate and curing area, 
z is the coordinate in the thickness direction, k* is 
the plate total thickness, g' is the transformed 

magnet ostrictive coupling modulus, fj is the 
magnetic field intensity, expressed in the following 
equation by Lee et al. (2004), 

H z (x,y,t) = k c I(x,y,t) with velocity feedback 

control / (x, y 9 t) = c(t) / in which k c is the coil 

dt 

constant, J(x 9 y 9 t) i s the c °ii current, c{t) is the 
control gain. 

Dynamic Equilibrium Differential Equations 

Without shear deformation effect, the time 
dependent of tangential displacement equations is 
assumed in the forms: u=u°(x 9 y 9 t) 

v = v°(x,y,^) and the transverse displacement 
equation of the middle-plane is assumed in the form: 
W = w(x 9 y 9 1) in which t is time. The dynamic 
equilibrium differential equations in terms of 
displacements including the magnetostrictive loads 
are expressed in the following matrix forms by Hong 
(2009): 

Ai 2 v4i 6 A 66 A 16 A 12 + A 66 A 26 0 0 0 

A 6 A 2 ^66 ^26 ^66 ^A 26 A 22 0 0 0 

B u +B l6 2 B l6 +B u +B 66 B 66 +B 26 B l6 +B 66 B u + B 66 +2B 26 B 26 +B 22 0 0 0 

f ay av av av av ay pw_ aV j* 

\ dx 2 dy 2 dx 2 5y 2 3x 2 5y 2 J 

\f\ | [10 0“ 

= < / 2 +p0 1 0 

Lo o -i 

Where f \ , f 2 , f 3 are the expressions of thermal loads (TV, M) , mechanical loads 9 p 29 q) and 
magnetostrictive loads ( TV, M ) , 
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^ fflV, 8N V 

1 dx dy 1 dx dy 

dN„ dN v dN ' dN v 

f 2 =—^ + ^ + p 2 +^^ + ^, 
dx dy ox dy 

J /"N /"N A ^ 

ox dy dy 

, afir, , 2 sm v , afir. 


0X 


<4 dy 


n 

(N X ,M X ) = (0i ]^ A + 0] 2 + 0i6^j)(^O5 Z T^^l)^ Z 


n 

(N y ,M y ) = J_2. (0 12 a x +g 22 « v +Q 26 a xy )(T 0 ,z-yT l )dz 


n 

(Vv > M xy ) = !_*• (016«, + 026«v + 06fi«xv )V >- ) dz 


n 

(iV x ,M,) = 0e J] // I (l,z , )&, 

2 

h*_ 

(N y ,M y ) = \ye 32 H z ( \,z 2 )dz, 

2 

(S,,M„) = £.<?*», (1. z ! )&, 

2 

(4 , 4 ,4/) = g. 4 a, z, z 2 )dz , (ij = l, 2, 6), 

2 

h*_ 

(P’H) = J_2. A(l,z)Jz. 

2 

in which yC> 0 is the density of ply, and /? 2 are 
the in-plane distributed forces, q is the applied 
pressure load. 


Dynamic Discretized Equations 

Without the shear deformation effect and under the 
following vibration of time sinusoidal displacement 
and temperature: 

u = [u°(x,y)\sm(co mn t), 

(8a) 

v = [v°(x,7)]sin(ffl m)i 0' 

(8b) 


w= w(x,y)sin(<y m „0, 

(8c) 

Ar = [r 0 (x,j) + -E7;(x,j)]sin(^). 

A 

(8d) 

where the term co mn is the natural frequency of 
plate, y is the frequency of applied heat flux. 

We apply the weighting coefficients of discretized 
equations (4) in the two-dimensional GDQ method to 
discrete the differential equations (7). And we use the 
following non-dimensional parameters under 
sinusoidal temperature 

(T 0 = 0, T x - T x sin(;r x! a)sin(;r y lb)). 

X — x! a, 

(9a) 

Y = ylb, 

(9b) 

U = u° I a, 

(9c) 

V = v° lb, 

_(9d) 

W = m*w/(cc x T l a 2 ). 

(9e) 

Thus, we can obtain the dynamic discretized matrix 
equations as in the paper by Hong (2009). 

Some Numerical Results and Discussions 

We like to consider the FGM plate consisting of two 
materials, of which the FGM material 1 is SUS304 
(Stainless Steel), the FGM material 2 is Si 3 N 4 

(Silicon Nitride). The temperature-dependent 
coefficients P Q9 P_ l9 P l9 P 2 and P 3 used to calculate 

material property terms, and E x E 2 V 1 V 2 p x 
p 2 CC X CC 2 and K x K 2 of these two typical 
constituent materials are listed in Table 1 by Shariyat 
(2008). The upper surface magnet ostrictive layer of 

the three-layer (0° ” / 0° /0°) laminates FGM plate 
under four sides simply supported is considered, the 
superscript of m denotes magnetostrictive layer. 
The elastic modules, material conductivity and 
specific heat of the Terfenol-D magnetostrictive 
material are used the same value as in the paper by 
Hong (2007). We use the grid points for the GDQ 
computation as in the following coordinates: 
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x. = 0.5[1 - cos(- — -7r)]a,i = 1,2,..., Nr 
1 W- 1 

(10a) 

y y =0.5[l-cos(^A*7 = l,2,...,M. 

M - 1 

(10b) 

The dynamic convergence results are obtained for 
center displacement amplitude w (a / 2, b / 2) 
without shear effects in the thermal vibration of 
sinusoidal temperature only (T 0 = 0, 7j = \§§ K , 
Pi = p 2 = q — 0 ) at time t = 6s, mode shape 
m — n — \ , with control gain value k c c(t ) = 0 , 
aspect ratio alb= 0.5, 1 and 2, side-to-thickness 
ratio alh = 100, 50, 20, 10 and 5, /z* =1.2 mm, 
h 3 = 1 mm, /z, = h 2 = 0.1 mm, /? ;/ = 1 , 

T = 300°K. Table 2 shows the w(a/2,b/2) 
(unit mm) in the grid point NxM = 17x17 
21x21 25x25 29x29 and 33x33 of GDQ 
method for the laminated Terfenol-D FGM plate at 
time t = 6 s. We find the NxM = 33x33 grid 
point have the good w(a/2,b/2) convergence 
result and use further in the GDQ analyses of time 

responses for displacement and stress for a I h = 5, 
20 and 100. 

We obtain the lowest frequency y of applied heat 
flux and vibration frequency CD n of Terfenol-D 
FGM plate ( m = n = 1 ), at time t = 0.001s, Is, 
2s,. ..and 9s, NxM= 33x33, h = 1.2 mm, 
/z 3 = 1 mm, h x = h 2 = 0.1 mm, alb = 1 , q = 0 , 

R n = 1, T = 653 K r T x = 100°^, as shown in 
Table 3. For the magnetic coil constant k c effect on 
the w (a / 2, b / 2) (unit mm) of vibration under the 
constant gain value c(/)= 1, time t= 6 s, alh =5 
without shear, we obtain the sketch of w(a / 2, b / 2) 
vs. k c as shown in Fig. 1. The suitable product 
value of k c and c(7) can be used to reduce the 
amplitude of w(a/2,b/2) into a smaller value 
near 0.00, there is an amplitude peak value at k c c{t ) 

= 7.16E08, we find k c c(t)= 1.56E09 for alh * =5 
thick plate can be chosen as the best suitable product 
values. 

Firstly, thermal sinusoidal vibration is investigated 
with time step equal to 0.1s, the suitable chosen 
product k c c{t ) of coil constant and controlled gain 


values versus time t for R n = 1 / thick plate 
alh =5 and thin plate a//z* =20,100 at 

T = 653 K , as shown in Table 4. 

Fig. 2 shows the w (a 1 2,6/2) (unit mm) versus 
time t of GDQ method for the laminated Terfenol-D 

FGM plate alh =5 , 20 and 100, respectively 
without shear effects. At time t= 0.001s, there is a 
great amplitude value of displacement with 

uncontrolled value ( k c c{t) = 0), w(a/2,b/2) = 

-0.156429 mm for thick plate alh =5 , 

w(a/2,b/2) = -4.31174 mm for thin plate 

alh =20, W (at 2,b/2) = -0.22975 mm for thin 

plate alh =100 . We find the amplitudes of 

w(a/2,b/2) with controlled k c c(t ) values are 

smaller than the amplitudes of w {a 1 2, b 1 2) with 
uncontrolled value (k c c{t}= 0), generally by using 

the GDQ method. We can use the suitable product 
values of controlled k c c(t ) to reduce the amplitude 

of w (a / 2, b / 2) into a smaller value near 0.00. 

Fig. 3 shows the time response of the dominated 
dimensional stress <j (unit GPa) at center position 

of upper surface Z = 0.5 h with respective to time 
for the laminated Terfenol-D FGM plate alh =5, 
20 and 100, respectively without shear effects 
a =0 and cr = 0 . We find the maximum 

xy xy 

response values of cr = -9.13E-04GPa for alh =5 
and 20, cj x = -9.12E-04GPa for alh =100 with 
controlled k c c{t ) case are almost equal to the 
response values of cr = -9.12E-04GPa with 

uncontrolled case ( k c c(t ) = 0), generally by using the 
GDQ method. 

Fig. 4 shows the compared non-dimensional 
W(X,b/ 2,6) versus X of alh* =5, 20 and 100 
without / with shear effects, for Terfenol-D FGM 
plate, at time t= 6s, m - n - \, NxM= 33x33 
for the case of without shear and 17x17 for the 
case of with shear, h = 1.2 mm, h 3 = 1 mm, h x = 
h 2 = 0.1 mm, alb = 1 , q = 0 , R n = 1 , 
T = 653 °K, T\ = 100°^, k c c(t)= 0. In the shear 

effect case, we use the YNS first-order shear 
deformation theory for the time dependent of 
displacement field, the value for shear correction 
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coefficients k a — kp — 516 used in the dynamic 
equilibrium differential equations by Hong (2012). 

For the thick plate at h = 5 , the values of 
displacement W(X,b/ 2,6) versus X without 
shear case are smaller than that with shear case. The 
maximum value W(X,b/2,6)= 0.0821655 occurs at 

X= 0.5 of the with shear case, and W(X,b/ 2,6) = 
-0.0325304 occurs at X = 0.691342 of the without 

shear case. For the thin plate a / h = 20, the values 
of displacement W(X,b/ 2,6) versus X without 

shear case are greater than that with shear case. The 
maximum value W(X,b/2, 6)= -0.00133522 occurs at 

X = 0.777785 of the without shear case, meanwhile 
W(X,b/ 2,6)= 0.00033884 occurs at 1= 0.5 of the 

with shear case. For the thin plate a / h =100, the 
values of displacement W(X,b/2,6) versus X 

without shear case are greater than that with shear 
case. The maximum value W(X,b/2,6) = 

0.0000138781 occurs at X = 0.0380602 of the without 
shear case, while W(X,b/ 2,6) = 0.000000541036 
occurs at X= 0.5 of the with shear case. The values 
of displacement W(X,b/ 2,6) are decreasing with 

a / h values increasing for square plate al b = 1 
with and without shear deformation effects. The 

deflection of center position versus al h at t= 6 
sec in Fig. 4d shows the deflection of thick plate 

a / h = 5 , and the case of deflection value with 
shear effect is much greater than the case of without 
shear effect. The deflections with and without shear 
effects are almost in the same values for the thin plate 
a/h* =100. 

Secondly, transient response is investigated with time 
step equal to 0.001s and use the fixed frequency 
y - 523.599/ S of applied heat flux. The suitable 
chosen product k c c(t ) of coil constant and 
controlled gain values versus time t for R n = 1 , 

thick plate a/h = 5 and thin plate a / h =100 
at T = 653 K, as shown in Table 5. We find the 
more thin plate is easier to control the displacement 
with less values of k c c(t) . 

Fig. 5 shows the transient value w (a / 2, b / 2) (unit 
mm) versus time t of GDQ method for the 

laminated Terfenol-D FGM plate a/h =5 and 100, 
respectively without shear effects. At time t= 0.001s, 
there is a great amplitude value of displacement with 


uncontrolled value ( k c c(t ) = 0), from w (a / 2, b / 2) 
= -0.156429 mm converges to small value for thick 
plate al h =5, from w {a 1 2,b/2) = -0.22975 mm 

converges to small value for thin plated / h =100. 
We find the amplitudes of w(a/2,b/2) with 

controlled k c c{t ) values are smaller than the 
amplitudes of w(a/2,b/2) with uncontrolled 
value ( k c c(t) = 0), generally by using the GDQ 
method. We can use the suitable product values of 
controlled k c c(t ) to reduce the amplitude of 
W {a / 2, b / 2) into a smaller value near 0.00. 

Fig. 6 shows the transient value of the dominated 
dimensional stress cj x (unit GPa) at center position 

of upper surface Z = 0.5 h versus time for the 
laminated Terfenol-D FGM plate a/h =5 and 100, 
respectively without shear effects a xy — 0 and 

<j xy — 0 . We find the transient values cj x oscillate 
between -9.12E-04GPa and 9.12E-04Gpa for 
a/h =5 and 100. The transient values cj x with 
controlled k c c{t ) case are almost equal to the 
values with uncontrolled case ( k c c{t ) = 0), generally 
by using the GDQ method. 

Conclusions 

The GDQ calculation provides a method to compute 
the controlled displacement and stress in the 

(0 m /0 /0 ) ply Terfenol-D FGM plate subjected 
to thermal vibration and transient response of 
sinusoidal temperature without shear deformation 
effect. The computation provides the following 
results, (a) The suitable controlled product values of 
coil constant and control gain k c c(t ) can be used to 

reduce the amplitude of center displacement 
W (a 1 2,b 1 2) into a smaller value near zero, (b) 
The amplitudes of stresses cj x of plates are almost 

in the same values under the two cases: with and 
without k c c{t ) values, (c) The values of 

non-dimensional displacement W(X,b/ 2,6) are 

decreasing with a/h values increasing for square 
plate a/b = l under the two cases: with and 
without shear deformation effects, (d) The deflection 

values of center position versus a/h at t= 6 sec 
are investigated, which found that the value with 
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shear effect is much greater than the case of without 

shear effect for the thick plate a / h = 5 , as well the 

deflections with and without shear effects are almost 

in the same values for the thin plate a / h =100. 
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fig. i w(a/2,b/2) vs. k c for a / h* = 5 , c{t)= l 
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FIG. 2a w(a/2,b/2) vs. t for a/h = 5 
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FIG. 2b w((l/2,b/2) vs. t for a/h = 20 


FIG. 3b C x vs. t for a/h =20 
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-^ uncontrolled (k c c(t) = 0) 
■ controlled (with^ c c(0 ) 


j t 


0123456789 


fig. 2c w(a/2,b/2) vs. t for a / h* = 100 
fig . 2 w(a/2,b/2 ) vs. t for a/h* =5,20andioo 

WITHOUT SHEAR 


FIG. 3c (J x vs. t for a I h* =100 

FIG. 3 (T x vs. t hr a/h* = 5 , 20 AND 100 WITHOUT 
SHEAR 
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W(X,b/2, 6) 



FIG. 4a W(X,b/ 2,6) vs. X for a I h* = 5 


0.0140 

0.0120 

0.0100 

0.0080 

0.0060 

0.0040 

0.0020 

0.0000 

- 0.0020 


I Meriinn 



FIG. 4d CENTER DEFLECTION vs. d/h* AT t = 6 sec 

FIG. 4 COMPARED W(X,b/2,6) vs.X ford/ h* =5, 
20 AND 100 




FIG. 4b W(X,b/ 2,6) vs. X for a! h* =20 


W(X,b/ 2, 6) 



fig. 4c W(X,b/ 2,6) vs - ^ for alh* =100 


FIG. 5a TRANSIENT VALUE OF w(a/2,b/2) vs. t FOR 

a! h* -5 



FIG. 5b TRANSIENT VALUE OF w(d / 2,b / 2) vs. 1 for 
d/h* =100 

FIG. 5 TRANSIENT VALUE OF w(d/2,b/2) vs. / for 
a! h* =5 AND 100 WITHOUT SHEAR 
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-1.5E-03 

-2.0E-03 


k ( c(l) = 0 


with k c c(t ) 



i i i i 

0 0.05 0.1 0.15 0.2 


FIG. 6a TRANSIENT VALUE OF <7 x vs. t FOR a / h* =5 


2.0E-03 

1.5E-03 

1.0E-03 

5.0E-04 

0.0E+00 

-5.0E-04 

-1.0E-03 

-1.5E-03 

-2.0E-03 


k c c(t) = 0 

with k c c(t) 



111 


.HI, 



II 


II 


0 0.05 0.1 0.15 '0.2 


FIG. 6b TRANSIENT VALUE OF (7 x vs. t FOR 

alh* =100 

FIG. 6 TRANSIENT VALUE OF (J x vs. t FOR Cl! h* =5 
AND 100 WITHOUT SHEAR 


TABLE 1 TEMPERATURE-DEPENDENT COEFFICIENT OF CONSTITUENT FGM MATERIALS 


Materials 

Pi 

P 0 

P - . 

Pi 

P 2 

p \ 

SU5304 

£1 (P a ) 

201.04E09 

0 

3.079E-04 

-6.534E-07 

0 


v , 

03262 

0 

-2.002E-04 

3.797E-07 

0 


p x {Kg / m 3 ) 

8166 

0 

0 

0 

0 


«, CK~ l ) 

1233E-06 

0 

8.086E-04 

0 

0 


k x (W lm°K) 

15.379 

0 

0 

0 

0 


C vl (J/Kg°K) 

496.56 

0 

-1.15 IE -03 

1.636E-06 

-5.863E-10 

Si 3 n 4 

P‘2 ( Pa ^ 

348.43E09 

0 

-3.70E-04 

2.16E-07 

-8.946E-11 


^2 

0.24 

0 

0 

0 

0 


p 2 (Kg/m 3 ) 

2370 

0 

0 

0 

0 


a 2 ( °K 

5. 8 723 E -06 

0 

9.095E-04 

0 

0 


k 2 ( W lm°K) 

13.723 

0 

0 

0 

0 


C v2 (J / Kg°K) 

555.11 

0 

1.016E-03 

2.92E-07 

-1.67E-10 
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TABLE 2 DYNAMIC CONVERGENCE OF TERFENOL-D FGM PLATE WITHOUT SHEAR 


n lh* 

GDQ method 


W {a / 2, b / 2) (unit mm) 


at n 

NxM 

alb = 0.5 

alb = 1.0 

al b = 2.0 

100 

17x17 

-0.634165E-04 

-0.905412E-04 

-0.175673E-03 


21x21 

-0.563660E-04 

-0.796810E-04 

-0.164627E-03 


25x25 

-0.518107E-04 

-0.735170E-04 

-0.16 029 IE -03 


29x29 

-0.494592E-04 

-0.709745E-04 

-0.158750E-03 


33x33 

-0.483320E-04 

-0.701375E-04 

-0.158280E-03 

50 

17x17 

-0.202713E-03 

-0. 28331 2E -03 

-0.585536E-03 


21x21 

-0.191526E-03 

-0. 27164 IE-03 

-0.579210E-03 


25x25 

-0.187638E-03 

-0.269502E-03 

-0.577589E-03 


29x29 

-0.186386E-03 

-0.269140E-03 

-0.577392E-03 


33x33 

-0.186190E-03 

-0.26951 IE-03 

-0.577425E-03 

20 

17x17 

-0. 99911 6E -03 

-0.131836E-02 

-0.205093E-02 


21x21 

-0.987938E-03 

-0.131330E-02 

-0.205050E-02 


25x25 

-0.981546E-03 

-0.131383E-02 

-0.205044E-02 


29x29 

-0.990436E-03 

-0.131042E-02 

-0.205094E-02 


33x33 

-0.990760E-03 

-0.130088E-02 

-0.20501 0E-02 

10 

17x17 

-0.241265E-02 

-0.283326E-02 

-0.293830E-02 


21x21 

-0.241384E-02 

-0.283509E-02 

-0.293858E-02 


25x25 

-0.24 115 IE-02 

-0.283366E-02 

-0.293846E-02 


29x29 

-0.241165E-02 

-0.283456E-02 

-0.293842E-02 


33x33 

-0.240542E-02 

-0. 28341 0E-02 

-0.29401 0E-02 

5 

17x17 

-0.340208E-02 

-0.343168E-02 

-0.232384E-02 


21x21 

-0.340130E-02 

-0. 34299 IE-02 

-0.232400E-02 


25x25 

-0.340583E-02 

-0.343044E-02 

-0.23238 7E-02 


29x29 

-0.340369E-02 

-0.343035E-02 

-0.232335E-02 


33x33 

-0.33 976 IE-02 

-0.343055E-02 

-0.232343E-02 



TABLE 3 

y AND CO n OF TERFENOL-D FGM PLATE WITHOUT SHEAR 



t 

= 0.001 s 


t= Is 


t= 2s 

al h 

r 

® ii 

r 

£0 n 

r 

® 11 

100 

523.599 

0.569829E-02 

1.57080 

0.141884E-01 

0.785398 

0.569830E-02 

20 

523.599 

0.567729E-02 

1.57080 

0.151387E-01 

0.785398 

0.567730E-02 

5 

523.599 

0.196036E-01 

1.57080 

0.136018E-01 

0.785400 

0.5381 13E-02 



t= 3s 


t = 4 s 


t=5 s 

al h 

r 


r 


r 


100 

0.523599 

0.569830E-02 

0.392699 

0.569830E-02 

0.314159 

0.569829E-02 

20 

0.523599 

0.567730E-02 

0.392699 

0.567730E-02 

0.314159 

0.567729E-02 

5 

0.523599 

0.5381 12E-02 

0.392701 

0.5381 13E-02 

0.314161 

0.5381 13E-02 

IT* ~ 


t= 6s 


t=7 s 


t= 8s 

al h 

r 

® 11 

r 

® 11 

r 

co n 

100 

0.261799 

0.569830E-02 

0.224400 

0.569830E-02 

0.196349 

0.569830E-02 

20 

0.261799 

0.567730E-02 

0.224400 

0.567730E-02 

0.196349 

0.567730E-02 

5 

0.261800 

0.5381 13E-02 

0.224401 

0.5381 12E-02 

0.196351 

0.5381 13E-02 


7 ^ 

100 0.174532 0.569830E-02 

20 0.174532 0.567730E-02 

5 0.174534 0.5381 13E-02 


21 



www.s e ipub .org/rame 


Researches and Applications in Mechanical Engineering (RAME) Volume 2 Issue 1, March 2013 


TABLE 4 k c c{t ) VERSUS t FOR R n =l WITHOUT SHEAR 


at h* 



k c c(t) 




t = 0.001- 0.1 

t = 02~ 0.4 

t = 0.5 ~0.9 

t = 1.0 t = 1.1- 1.4 

f=1.5 

5 

-10 9 

-3xl0 6 

-10 7 

ON 

o 

1 

C4 

o 

1 

-10 7 


t = 1.6~ 3.3 

? = 3.4~ 3.9 

NO 

I 

Q 

II 

II 

X 

NO 

o 



-10 9 

-2xl0 9 

10 9 

1.56xl0 9 



t= 0.001 

? = 0.1 t = 0.2 

t = 03-0.9 

t= 1.0 7=1.1 

7=12- 1.3 


-10 9 

ON 

o 

1 

NO 

o 

-10 5 

ON 

o 

1 

<N 

o 

1 

-10 7 

on — 

t = 1 4~ 3.2 

? = 3.3~ 3.5 

t = 3.6 

t = 37-5.2 t 

= 53-5.4 

zU 

10 9 

10 5 

10 2 

-10 9 

10 5 


t = 5.5- 7.5 

t = 7. 6~ 9.0 





-10 7 

10 9 





t = 0.001 

? = 0.1 ~ 1.5 

t = 1.6 ~1.9 

7 = 2.0 t = 2.1~ 2.3 

t = 2.4 

100 

-10 9 

10 2 

-10 9 

-10 2 10 7 

-10 9 


t = 2.5- 9.0 
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